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D ■ ABSTRACT 
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I review galaxy number count models and present ncmod, an extendable and 

^ ■ general purpose model for comparing and interpreting the results of field galaxy 

P!!^ ■ survey data. I develop techniques and software for converting the results of a 

O ' survey done in one filter into another filter, for direct comparison with other 

^ ■ surveys. Comparison of the data from surveys which differ greatly in wavelength 

OO . coverage or sensitivity is of necessity model-dependent, but comparison between 

. similar surveys can be done in a relatively model-independent way. I extrapolate 

P_|. existing number counts into the ultraviolet and thermal infrared. The model 

Q . is used to predict the results of future space missions, including STIS and 

U NICMOS on HST, ISO, SIRTF and NGST. 
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^ ■ 1. Introduction 

Counting the surface density of galaxies on the sky as a function of apparent magnitude 
is one of the classic cosmological tests. Both the power and the difficulties of extragalactic 
galaxy surveys were recognized by Hubble in the 1930s, when he attempted to constrain 
the cosmological geometry through galaxy counts, but was limited by the difficulty of doing 
accurate galaxy photometry with the poor quality data he had available (Hubble 1934; 
Hubble 1936), and the inaccuracies of his K-corrections. In the 1960s, in a classic review of 
observational tests of world models possible with the Mount Palomar 200-inch telescope, 
Sandage (1961a) included galaxy counts as one of four fundamental tests of observational 
cosmology. He noted, however, that the number- redshift relation n{z) is more sensitive to 
the value of go than the number-magnitude relation n{m). He noted further that at the 
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limit of the 200-inch telescope, which he took to be = 22, the difference in the number 
counts measured in an open universe with = and a closed universe with qq = 0.5 is just 
AniR — 0.28, which was comparable to the photometric and systematic errors of galaxy 
photometry of the time. He claimed that it would be impossible to find the correct world 
model from galaxy counts, a statement that is arguably true even today when we have 
counts Imag deeper than the limit he considered. To the reasons he listed, systematic and 
statistical photometric errors, aperture and surface brightness corrections, and uncertainties 
in the idealized models, must now be added the unknown details of galaxy formation and 
evolution. However, some of the more extreme cosmological models which he considered, 
including those with q'o > 1 or the steady-state model with — —1, can be addressed with 
existing number count data. 

Although also discussed by Sandage (1961b), the effects of luminosity evolution of 
galaxies on the n{m) relation, were not studied in detail until the 1970s. In a series of papers 
culminating in a review of her galaxy count model (Tinsley 1977; and references therein) 
developed the technique of population synthesis, in which she took an initial mass function 
(IMF) of stars, and observed broadband stellar colors to predict the color and luminosity 
evolution of stellar populations. To construct her model, she made several assumptions 
which would dominate the study of galaxy counts for the next two decades. She assumed 
that all galaxies formed at the same time, and that the time-scale of star-formation was the 
factor which determined the relation between morphological Hubble type and SED. In this 
picture, elliptical galaxies formed their stars in a short initial burst, while spiral galaxies 
had longer timescales of star-formation. She assumed that there was conservation of galaxy 
number, and merging was not an important factor. 

The population synthesis approach dominated studies of galaxy evolution in the 
1980s, when workers used stellar spectra instead of broadband colors to construct galaxy 
templates. In a series of papers Bruzual (1983a) extended the technique to the ultraviolet 

(UV) using new observations of stars with the lUE telescope. At about the same time, 
Arimoto & Yoshii (1986), 1987 and Guiderdoni & Rocca-Volmerange (1987) included the 
effects of chemical evolution of galaxies, tracing the change in metallicity as the population 
of stars within the galaxy evolves. These results were included in galaxy count models by 
Yoshii & Takahara (1988) and Guiderdoni & Rocca-Volmerange (1990). 

1.1. Galaxy Count Models 

The galaxy count model of Yoshii & Takahara (1988) serves as the basis for this paper, 
because many of the ingredients for the model have now become standard. They began 
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with a measured local luminosity function. They assumed that the characteristic magnitude 
(M* in the Schechter 1976 parameterization) and faint end slope (a) of the himinosity 
function was independent of Hubble type in the bj filter, and set the distribution of types 
(and thus 0* as a function of type) to match those observed in a local survey. They then 
converted that luminosity function to other filters using type-dependent local galaxy colors 
from their population synthesis model. This has the effect that while a is the same for 
each type of galaxy in each filter, M* is a function of type in filters other than bj, and the 
summed luminosity function in other filters has an overall a that is different from that of 
the individual types. Examples of this effect are given below. Other galaxy count modellers 
(Guiderdoni & Rocca-Volmerange 1990) used type-dependent luminosity functions (e.g. 
King & EUis 1985), even though these were of necessity less accurately determined. 

The population synthesis models give the SED of each type of galaxy as a function 
of time since the galaxy formed. (Following Tinsley 1977, the morphological galaxy types 
are each equated with a star-formation history). The chosen cosmology give the relation 
between redshift, distance modulus and K-corrections, and the time-redshift relation is 
applied to the luminosity evolution model so as to include evolutionary corrections. The 
model is then integrated over redshift to give the n(m) relation, or integrated over other 
parameters to give n{z), and color distributions. 

Recently a new type of model has appeared. The semi-analytic approach, motivated 
by cosmological theory, was proposed by White & Frenk (1991) and Cole (1991), and 
developed further by Lacey et al. (1993), Kauffman, White & Guiderdoni (1993), Kauffman, 
Guiderdoni & White (1994) and Cole et al. (1994). This type of theory begins with the 
physics of the big bang, cold dark matter, and hierarchical structure formation, and traces 
the formation of galaxies ab initio using a semi-analytic approach constrained by numerical 
simulations. It then predict observables, such as the number counts, luminosity functions, 
colors, redshift distributions, and the TuUy-Fisher relation. Baugh, Cole & Frenk (1996) 
have extended this model to study the distribution of morphological Hubble types, as 
measured by the bulge-to-disk ratio. The fundamental difference between this type of 
model and the traditional number count models is the direction in time. Traditional models 
take the local population and attempt to extend it backwards in time, and determine 
what it evolved from. The semi-analytic approach takes the initial conditions set by 
the big bang and attempts to trace their evolution into the local population. While the 
semi-analytic theory is the most physically motivated, and may eventually provide the 
greatest understanding of galaxy formation and evolution, there is still a place for the 
traditional models. Traditional models are still the best way of comparing different data 
sets, and this is the approach I take in this paper. 
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1.2. Galaxy Count Data 

A considerable amount of observational galaxy count work was done in the 1980s using 
photographic plates, both on Schmidt telescopes to get bright counts, and on 4m class 
telescopes to get faint counts. The modern era of galaxy count data, using linear detectors 
(CCDs) to count galaxies to very faint levels began with the work of Tyson (1988) and 
Lilly, Cowie & Gardner (1991). Those workers estabhshed the counts in several optical 
filters to S ~ 28, while a series of papers by Metcalfe et al. (1991) estabhshed CCD counts 
at intermediate magnitudes. These data, when compared to the modeling work described 
above, led to what has become known as the faint blue galaxy problem. The data showed an 
excess in the number of galaxies over the no-evolution model predictions. While the initial 
interpretation was that the surveys were detecting the ultraviolet light from star-forming 
galaxies in an open cosmology, further investigation, including redshift surveys (Broadhurst, 
EUis & Shanks 1988; CoUess et al. 1990; Cowie, Songaila & Hu 1991) showed that this 
excess was not at high redshift. The excess was not seen in the near-infrared band 
number counts (Gardner, Cowie & Wainscoat 1993), and was made up of blue galaxies. 

The spatial resolution of images taken by the Hubble Space Telescope (HST) has 
allowed the morphological study of faint galaxies. The HST Medium Deep Survey 
key project (Griffiths et al. 1994; hereafter MDS) has proven to be a great success 
in characterizing the optical morphology of galaxies at intermediate magnitudes. In 
conjunction with the fainter images of the Hubble Deep Field (WiUiams et al. 1996; 
hereafter HDF), the MDS has revealed an excess of galaxies with unusual or disturbed 
morphologies, over what is predicted by the number count models (Glazebrook et al. 1995; 
Driver, Windhorst & Griffiths 1995; Abraham et al. 1996). These irregular galaxies likely 
make up the population of excess faint blue galaxies. 

Bright 6 J— band galaxy counts, measured on scanned photographic plates covering 
4300 square degrees, show a steeper slope than is predicted by the models (Maddox et al. 
1991). This steep slope was initially believed to be due to rapid evolution at low redshift, 
but recently the data has been called into question. Systematics in the photometry, 
either due to the rapid plate scanning process, (Metcalfe, Pong & Shanks 1995), or to the 
distribution of central surface brightness in the galaxies, (McGaugh 1994), could steepen 
the slope of the number counts artificially. This has caused considerable debate about the 
magnitude of the faint blue galaxy problem. When the models are normalized at the bright 
end of the local steep slope, B ~ 15, the excess in the counts is a factor of 2 to 4 over the 
model predictions at i? ~ 22. The results of redshift surveys of galaxies selected at this 
magnitude level, however, show little departure from no-evolution predictions of the shape 
of the number-redshift relation (CoUess et al. 1990), although their absolute numbers, or 
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normalization showed the excess seen in the number counts. If the models are normalized 
at the faint end of the local steep slope, B ~ 18, and passive evolution is included in the 
model predictions, then the faint blue galaxy excess does not appear until B > 25 for the 
closed universe case, which is fainter than any of the spectroscopic redshift surveys. Recent 
bright number counts, obtained with CCD imaging of 10 square degrees, do not confirm the 
steep slope, and are consistent with the higher normalization (Gardner et al. 1996). Larger 
CCD surveys are needed to confirm this result. 

There have been several attempts to modify the traditional number count models 
to explain the faint blue galaxy problem by introducing one or more additional free 
parameters into the model. Rocca-Volmerange & Guiderdoni (1990) and Broadhurst, 
EUis & Glazebrook (1992) relaxed the constraint of conservation of galaxy number, and 
introduced large amounts of merging. A disappearing dwarf galaxy model (Dekel & Silk 
1986; Cowie et al. 1991; Babul & Rees 1992) has gained support from deep redshift 
surveys (Lilly et al. 1995; Ellis et al. 1996) which see a steepening of the faint end slope 
of the luminosity function to higher redshift, particularly for the bluer galaxies. High 
normalization models have much less of a faint blue galaxy excess to explain, and succeed in 
fitting most of the observations with modifications to the luminosity evolution models. One 
recent attempt used a steep slope for the luminosity function of late-type spiral galaxies, 
which has an effect on the counts similar to the disappearing dwarf models (Metcalfe et al. 
1996). Cowie et al. (1996) recently proposed a downsizing model, in which less luminous 
galaxies evolve rapidly at low to intermediate redshift, while bright galaxies evolve rapidly 
at redshift z > 1. 

These models have usually attempted to fit all of the available field galaxy survey 
data, and do so by introducing additional free parameters into the models. Fitting the 
data, however, is often only done to within a factor of ~ 2. In this paper I take a more 
modest approach. Rather than trying to construct a single model which fits all of the 
data, I construct an extendable or general purpose model. The model can be tuned to 
fit one observation, and then used to predict an observation in another filter. This is 
essentially a translation program, translating observations between filters. Number counts 
and luminosity functions measured in different surveys can be directly compared. When 
the filters are similar, such as the ground-based U band filter and the WFPC2 L^soo filter, 
this can be done in a relatively model independent way, as a model which fits the [/—band 
number counts is likely to also fit the C/300 counts. When the filters have very different 
central wavelengths, such as B and K, or when the observations are at very different 
magnitudes, then the comparison will depend more strongly on the model parameters 
chosen. I also will use this program to predict the results of future ground and space 
mission surveys, by extrapolating the most relevant existing data into the filters and depths 
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that will be achieved. 

In section 2 I will outline the equations and parameters that make up the model. In 
section 3 I will use the model to compare existing survey data, and to predict the results of 
future surveys. Section 4 is a summary. 

2. The Number Count Model 

The components which go into a number count model include the spectral energy 
distributions (SEDs) of the galaxies, as a function of galaxy type, either local SEDs to 
produce a no-evolution model, or SEDs that are a function of time to produce an evolution 
model. Next, a description of the local population of galaxies, including the distribution of 
galaxy types, and the luminosity functions for each type. Finally, there is the cosmological 
model, or the choice of qo, Aq and Hq. The model is then integrated. 

2.1. Spectral Energy Distributions 

The number count model begins with the spectral energy distribution (SED) of 
galaxies. The main method I use to make this model extendable is in the number and 
distribution of spectral types of galaxies. While previous models have used 3 to 5 spectral 
types, and equated them to morphological Hubble types, (e.g. E/SO, Sa, Sb, Sc, Irr, 
Yoshii & Takahara 1988), in this model I allow for any number of types, and include 
the ability to interpolate between the types to match an observed color distribution. 
The distribution can be determined by one set of observations, and then the model can 
be used to predict another set of observations. Both theoretical and observed SEDs 
are available in the literature, the theoretical SEDs coming from population synthesis 
models described above. Bruzual & Chariot (1996, in preparation) have made their galaxy 
isochrone synthesis spectral evolution library (hereafter GISSEL96) model available at 
ftp://gemini.tuc.noao.edu/pub/charlot/bc96, and in Leitherer et al. (1996). This model is 
ideal for the purposes of this paper because it takes the star-formation history of the galaxy 
and returns the SED as a function of time. Any star-formation history is possible, and so 
it is possible to construct a large number of galaxy spectral types. They have included the 
effects of differing metallicities, and find that galaxies have an age-metallicity degeneracy 
in their colors. 

Kinney et al. (1996) have taken an observational approach to determining the SEDs 
of galaxies. Using lUE pointed and archival data, and ground-based optical spectra on 
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the same galaxies through a similar aperture, they have obtained spectra from 1200A to 
lOOOOA for 30 quiescent and 36 star- bursting galaxies, with a range of morphological type. 
By averaging together the spectra of each type, they present a library of spectra for 5 
quiescent types, the bulge component, and 6 starbursting types. The starbursting types 
are differentiated by their intrinsic color excess E{B — V), caused by dust obscuration, 
which is the dominant parameter that determines the appearance of the starburst SED 
(Calzetti, Kinney & Storchi- Bergman 1994). The SEDs have been made available at 
f t p : / /f t p . stsci . edu / catalogs /spectra . 

Most of the Kinney ct al. (1996) SEDs extend to Ifim wavelength, but to match the 
extensive ground-based observations done in the K band at 2.2/im (e.g. Gardner et al. 
1993), to match the work done in the far-IR with IRAS (e.g. Hacking & Soifer 1991), 
and to predict the results of planned space missions in the near-, mid-, and far-IR, it is 
necessary to extend the SEDs to other wavelengths. Schmidt et al. (1997) have done this, 
by using previously published data available through the NASA Extragalactic Database 
(NED) to compile the SEDs of the Kinney et al. (1996) galaxy sample over a large range 
in wavelength. While they attempted to correct for the very different apertures used in 
different observations, the processes affecting the SED of galaxies vary widely at different 
wavelengths, and make it difficult to use their SEDs for the determination of K-corrections. 
Devereux & Hameed (1997) have shown that there is not a simple relation between galaxy 
type and far-IR SED, and each type of galaxy shows a large spread in its ratio of Lfir/ Lopt- 
My model allows the use of any number of SEDs, and allows interpolation between different 
SEDs, so it is possible to use a range of far-IR fluxes for each galaxy template. For the 
remainder of this paper, I will use SEDs from the GISSEL96 models (which include only 
stellar emission, and thus are not readily applicable in the Far-IR), however, my model can 
be easily extended to include non-stellar emission. This extension will be considered in a 
future paper. 

2.2. Luminosity Functions 

There now exist measurements of the luminosity function of field galaxies in the B 
(Lovcday ct al. 1992; Marzkc, Huchra & Gcllcr 1994a), R (Lin et al. 1996), and K (Gardner 
ct al. 1997) bands, and in the far-IR (Saunders ct al. 1990; Isobe & Feigelson 1991). In 
the optical and ncar-lR, a Schcchtcr (1976) parameterization is often used, while in other 
bands, a Gaussian or power-law/Gaussian combination is used. The Schechter (1976) 
parameterization is 

(t){L)dL = (t)*{L/L*Yexp{-L/L*)d{L/L*), (1) 
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or equivalently 

iP{M)dM = Mi^0*ea;p{-^^^(a + 1)(M - M*) - exp[-^^^(M - M*)]}dM, (2) 
2.5 2.5 2.5 

where exp{x) — e^, ln{x) is the natural logarithm, a is a constant index, 0* is a 
constant coefficient which has the dimension of the number density of galaxies, L* is the 
characteristic luminosity, and M* is the corresponding characteristic magnitude. Generally, 
the measurement of a luminosity function is also done in a non-parametric manner 
(Efstathiou, Ellis & Peterson 1988) and a discrete luminosity function is returned. 

There have been several determinations of the local luminosity function of galaxies 
measured as a function of type (Loveday ct al. 1992; Marzke et al. 1994b), color (Metcalfe 
et al. 1991; Shanks 1990; Marzke & da Costa 1997), or emission line properties (Lin et al. 
1996), as well as studies of the evolution of the LF with redshift (Lilly et al. 1995; Ellis 
et al. 1996; Cowie et al. 1996). The consensus of these studies is that red, early- type, 
and non-emission-line galaxies tend to have a shallower slope, (i.e. a > —1.0), than the 
blue, late-type, and emission-line galaxies {a < —1.0). While there is a consensus on the 
trend, however, there is not a consensus on the actual measurement of a, which is highly 
correlated with M* in any case. I have included the capability within my model to use a 
different LF for each spectral type. 



2.3. Cosmological Parameters 

Number count models are constructed in co-moving coordinates, as that is how 
luminosity functions are formulated. Following Yoshii & Takahara (1988), the relation 
between apparent magnitude and absolute magnitude in a filter F is: 

mp = Mf + Ef{z) + 5log{DL/l0pc), (3) 

where log is the base 10 logarithm, and 

- /„~A.te(0))f(A')dA' • 

where F{\) is the filter throughput. In the no-evolution case, Ep{z) is evaluated at taiO), 
and is called the K-correction. 

In the Friedman-Robertson- Walker model, the luminosity distance -Dl is defined as 
(Carroll, Press & Turner 1992; Fukugita et al. 1990): 

= jn^'^^^il^'^l'^' /'[(I + + 2(go + Ao);^') - z'{2 + z')Xo]-'/'dz'}, (5) 

-"0 "fc ' "'O 
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where Hq, q^, and c are the Hubble constant, the deceleration parameter, and the speed 
of light, respectively, and Aq is the normahzed cosmological constant, Aq = Aqc^/3Hq. 
is the curvature term, = 1 — 2go — 3Ao, and sinn is defined as sink for models with 
fifc > (open universe) and sin for models with < (closed universe). In a fiat universe, 
flk = 0, and the sinn and fi^s disappear from the equation, leaving only the integral. The 
comoving volume V{z) and the cosmological time t{z) differentiated with respect to z are 
given by 

dV AncDl 

(6) 



and 



dz Ho{l + zy[2{qo + Ao)(l + z)^ + Aq + l^fc(l + z^y/^ ' 
dt -1 



(7) 



dz Hoil + z){[2{qo + Xo)z + 1 - Ao](l + z)"^ + Aoj^^ ' 
respectively. The age tdz) of a galaxy formed at redshift zp is obtained by integrating 
dt/dz from zp to z. Hq cancels out in this model in every factor except tciz), and thus has 
no effect on no-evolution models. In models with evolution, a change in Hq is degenerate 
with a change in zp, as long as the current age of the galaxies is less than the Hubble time. 



2.4. Constructing the Model 

The galaxy number count data are obtained by counting all the images of galaxies 
in a finite area of the sky. If n{mx, z)dm\dz is the number of galaxies between m\ and 
mx + dm\ and between z and z + dz, then for (0 < z < zp), 

n{mx, z)dm\dz = — '^^ip{m\, z)dm\dz, (8) 
47r dz 

where uj is the angular area in units of steradians over which the galaxies are counted, 
and ip{mx, z) can be obtained by solving equation ^ for Mx (for each galaxy type) and 
putting the result into equation ||. I will use 1 deg"^ for our area, which corresponds to 
UJ = 3.05 X 10~^sr. Integrating n{mx, z)dmxdz with respect to z, gives the differential 
number count 

n{mx)dmx= / n{mx, z)dmxdz, (9) 
Jo 

which is summed over the distribution of galaxy types. 



2.5. Model dependence 



There are a considerable number of free parameters in this model, although some can 
be determined from local survey observations. In this section I consider the effects on the 
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model of varying the free parameters one at a time, while holding the other parameters 
fixed. The fiducial model I will use for comparison is a model of the number counts in the 
WFPC2 Jgu filter, for a model with go = 0.5, Hq = 50 km sec'^ Mpc-\ and Aq = 0.0, 
including the effects of evolution. I consider the effects of varying the galaxy mix below, 
so in this section I will use a galaxy mix consisting of 4 types with a single redshift of 
formation, and exponential rates of star formation with e-folding times of 1 Gyr, 4 Gyr, 7 
Gyr, and constant star formation, corresponding roughly to E/SO, Sab, Sbc and Scds. Half 
of the E/SO types have metallicities of 2.5 times solar, the other half have solar metallicity. 
The Spiral galaxies have solar, 2/5, and 1/5 solar metallicity respectively. The galaxies were 
formed at Zform = 15, and are in a distribution of 32%, 28%, 29%, and 5%, respectively. 
I include one additional starforming dwarf type, with constant starformation, 1/5 solar 
metallicity, and an age of 1 Gyr at every redshift. This last type does not passively evolve, 
and represents a population of starforming galaxies that is continuously refreshed at all 
redshifts (Gronwall & Koo 1995; Pozzetti, Bruzual & Zamorani 1996). For this fiducial 
model, I have used the type-independent i^'— band luminosity function as measured by 
Gardner et al. (1997), (M^ = -24.62, a = -0.91, and 0* = 2.08 x 10"^), converted to the 
type-dependent /814 band as described above. The models are plotted in Figure ^ along 
with a compilation of the /—band data. The galaxy types are plotted separately to show 
the contribution of each type to the total. 

In the past, number count model predictions have often been plotted as a function of 
cosmology, and with and without the effects of evolution. I do this in Figure |l|. There are 
several aspects of this fiducial model to note. First, the no-evolution models underpredict 
the number counts at the faint end in all of the cosmologies plotted. Second, evolution 
causes the model to shoot up rapidly at /814 ~ 20 for the go = 0.5 case, and at J814 ^ 22 
for the other two cosmologies. This effect, although quickly swamped by volume effects in 
the go = 0.5 model, results in an overprediction of the observed counts, and is due to the 
rest-frame ultraviolet flux of the massive stars involved in the initial burst of star-formation 
of the galaxies. This overprediction was absent from many early models. Some were based 
upon the tabulated Bruzual (1983b) models which only extended to 2 = 2, and thus did 
not include the effects of this initial burst of star-formation. Others used a parametric 
form for star- formation (i.e. Tinsley 1977), which also under-estimated the initial burst of 
star-formation. Wang (1991) proposed that extinction of the UV light by dust produced 
by the star-formation (which is not included in the GISSEL96 models) will eliminate this 
overprediction of the models (Gronwall & Koo 1995; Babul & Ferguson 1996; Roche et al. 
1996). I consider this effect below. In addition, Pozzetti et al. (1996) used a Scalo (1986) 
initial mass function, rather than a Salpeter (1955) IMF to reduce the production of UV 
light. 



- 11 - 



The models plotted in Figure ^eacli have Zform = 15, and I now consider the effects on 
the models of varying the redshift of galaxy formation. In Figure I plot the results of the 
models run with Zform = 5 and -Z/orm = 2.5. The redshift of formation causes a turnover in 
the number counts. At fainter magnitudes, the surveys are no longer reaching to higher 
redshift, and therefore the number counts are no longer going up due to volume effects. The 
limiting slope of the number counts is set by the faint end slope of the luminosity function. 
This has led some workers to propose that there is a component of the local luminosity 
function with a steep slope, which comes to dominate at the faint end (Driver et al. 1994a). 
There is evidence that dwarf elliptical galaxies have a steep slope in the luminosity function 
measured in some clusters (Sandage, Binggelli & Tammann 1985; Driver et al. 1994b). In 
the field there have been attempts to detect this population, (Marzke et al. 1994a; Loveday 
1997), but it is not yet clear whether it is a universal phenomenon, or a feature of overdense 
regions. 

2.6. Internal Absorption by Dust 

In this section I consider the effects on the models of absorption by dust internal to the 
galaxies. First investigated by Wang (1991), it was recognized that the population synthesis 
models predicted too much UV flux in the absence of dust (Gronwall & Koo 1995; Roche 
et al. 1996; Babul & Ferguson 1996). Following Bruzual, Magris & Calvet (1988), Wang 
(1991) modeled the dust internal to a galaxy as an absorbing layer, symmetric around the 
midplane of the galaxy, whose thickness is a fraction ( of the total thickness of the stellar 
disk. The observed luminosity is 

Lobserved = L^l -()(! + e-^)/2 + ((1 - e-^)/r], (10) 

where r is the optical depth of the absorbing layer. I take ( = 0.25, and assume that the 
optical depth r oc L^, where P = 0.5. I will assume that the extinction is a power law in 
wavelength, oc A~", with n = 2. I set the present day extinction of galaxies to be r = 0.2 
for L* galaxies at 4500A, and scale all other galaxies according to the above equations. 
Using the assumption that the measured luminosity function has not been corrected for the 
effects of internal absorption by dust, the rest-frame correction is subtracted. 

Figure ^ is a plot of the model prediction for the /814 number counts, including the 
effects of dust. In comparison to the models plotted in Figure |1|, the models no longer 
overpredict the counts for the Aq = models, as the UV flux from the first burst of star 
formation is no longer seen at bright magnitudes. Further evidence for the role of dust in 
galaxy formation comes from long-slit spectroscopy and narrow-band imaging searches for 
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line emission from starbursts at high redshift (Thompson & Djorgovski 1995), and from 
models of damped Lya systems(Pei & Fall 1995). 

2.7. Merging, or Galaxies in Pieces 

Following suggestions by Koo (1990) and Guiderdoni & Rocca-Volmerange (1990), 
Rocca-Volmerange & Guiderdoni (1990) proposed that the galaxies in the present-day 
universe were formed in pieces which subsequently merged together. In their model, this 
number evolution takes the form 0* oc (1 + z)^ in the Schechter parameterization, and in 
order to conserve the luminosity density, L* oc (1 + 2;)"''; here, ?7 is a free parameter. This has 
the effect of dividing the flux from a single galaxy at high redshift between several galaxies. 
While this is usually considered to be due to the physical merging of galaxies, it also would 
be the result of an observational artifact introduced into the data by over-enthusiastic 
deblending in the object detection routine. CoUey et al. (1996) recently proposed that this 
observational effect is operating in the catalogs of the HDF, and is caused by the fact that 
galaxies tend to be clumpy and irregular in the rest-frame UV. Broadhurst et al. (1992) 
proposed a slightly more complicated function of redshift, 0* oc exp{—Q / [3{{1 + z)~^ — 1)), 
where Q defines the merger rate, and /5 is a function of the look-back time. This function 
avoids the unreasonably high merger rate at high redshift of the exponential form, and has 
the intuitive advantage that Q is approximately the number of pieces at z ~ 1 that merge 
to form a present-day galaxy. Glazebrook et al. (1994) used a simplified version of this 
function, with the merger rate oc 1 + Qz. Guiderdoni & Rocca-Volmerange (1991) give a 
detailed analysis of self-similar merging scenarios with the assumption of conservation of 
total comoving luminosity density, which leads to: 

L\z)<p\z)T{a{z) + 2) = L*(O)0*(O)r(a(O) + 2), (11) 

where r(a + 2) is the Gamma function that results from integrating the Schechter luminosity 
function. They alternatively hold each of the Schechter parameters constant while varying 
the other two. In this paper I will only consider merging models with a held constant, but 
the model is easily modified for the case when a is a function of redshift. Both Broadhurst 
et al. (1992) and Carlberg (1992) include in their models the effects of merging on the 
luminosity evolution of the galaxies, but I will not consider that here. 

Figure ^ is a plot of the model predictions for the fiducial model with Rocca-Volmerange 
& Guiderdoni (1990) merging, and with Broadhurst et al. (1992) merging. The two models 
give very similar results. In general, the decrease in luminosity of individual galaxies 
moves to fainter magnitudes the sharp upturn in the counts due to the UV flux from the 
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initial burst of star-formation. The increase in number density causes the counts at fainter 
magnitudes to rise. The go = 0.5, Aq = 0.0 cosmology is no longer ruled out by the high 
counts at the faint end. 



2.8. Excess Dwarf Models 

Number count models with an excess of dwarf galaxies have also been proposed as a 
way of reconciling the high counts with a go = 0.5, Ao = 0.0 cosmology. They have also 
been used to explain the number-redshift relation (Cowie et al. 1991), and the observed 
excess of irregular galaxies seen in the HST Medium Deep Survey (Glazebrook et al. 1995; 
Driver et al. 1995), and in the Hubble Deep Field (Abraham et al. 1996). The dwarf 
galaxies in these models are actively star-forming, a requirement introduced to fit the 
increase in emission line strength seen in the spectroscopic redshift surveys. I include this 
type of galaxy population in my model, and plot the results in Figure ^ I give the constant 
star- formation galaxy type a steep slope a = —1.8 in the luminosity function, balanced by 
a fainter = —23.12, and 0* = 7.5 x 10^'^. The dwarf galaxies have a negligible effect 
at bright magnitudes because of the faint M*, but dominate the number counts at faint 
levels due to the steep faint end slope. This model is somewhat ad hoc, as are many of 
the modifications to the standard models. It allows any arbitrary faint-end slope for the 
number counts, with the only stipulation being that the local dwarf population be too faint 
to be observed in the local redshift surveys used to determine the luminosity function. The 
faint end slope of the number counts is not as steep as this model predicts, so I also plot a 
model where the dwarf galaxies are given a faint end slope of a = —1.5. 

Figure |^ is a plot of the model predictions for dwarf galaxies plus dust, and for 
Rocca-Volmerange & Guiderdoni (1990) merging plus dust. This figure illustrates the 
interaction between the free parameters of the model. 

3. Applications of the Model 

In the last section I reviewed the ingredients of a traditional galaxy count model. By 
including all of the extra free parameters that have been invoked to fit the data, I have 
made it possible to fit observations of galaxy counts in a variety of different ways. I have 
shown that contrary to some early hopes, it is not realistic to expect galaxy counts alone 
to constrain the cosmological geometry, or even to strongly constrain the form of galaxy 
evolution. While some values of the free parameters are physically more realistic than 
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others, nonetheless there is too much degeneracy in the possible models. Surveys, however, 
contain much more data than just the number of galaxies. Most photometric survey fields 
have been observed with more than one filter, giving color information on the galaxies. 
Galaxy evolution is measured by studying the change with redshift of the general luminosity 
function, and spectroscopic redshift surveys provide the greatest leverage for these studies. 
Models like the one 1 have presented here can be used to directly interpret survey results. 

My purpose in this paper is not to construct a single model which fits all the data, 
however. The power of galaxy count models like this one is in comparing different 
observations with each other. The simplest application is a consistency check: e. g. is 
the local 6j— band luminosity function, as measured by Loveday et al. (1992), consistent 
with the faint S— band number counts, as measured by Tyson (1988). Answering this 
question in the negative led to the faint blue galaxy problem, and a modification of the 
models. In this section I will extend this use of the model to comparing observations made 
in different filters. I will use the model to convert measurements of the galaxy counts, and 
of the luminosity functions from one filter to another. This will serve two purposes. First, 
observations made in different filters can be directly compared with one another. Second, it 
will be possible to make predictions for the results of future observations using new filters, 
especially for future space missions that will address wavelength regimes for which there is 
very little existing data. 



3.1. Conversion between filters 

Comparing the results of a survey conducted in one filter to a survey conducted in a 
different filter can be troublesome. While previously converting number counts between 
similar filters has been done with a constant offset in magnitude (Metcalfe et al. 1991), 
this doesn't work for comparing observations between filters that are different enough 
that the slope will change. It is possible, however, to use a number count model to 
translate or interpolate measured counts from one filter to another. This is of necessity 
model-dependent, but when translating counts between very similar filters, such as from 
the Johnson Uj filter to the WFPC2 U300 filter, or from the ground-based K filter to 
the NICMOS Hi q filter, the choice of model parameters makes little difference. When 
translating counts between filters with widely separated central wavelengths, or when 
extrapolating ground-based counts into the UV or thermal IR, the resulting prediction 
depends more strongly on the model used. 

To transform a measurement of the number counts from one filter to another, 1 first 
find a model which matches the counts in that filter reasonably well. When interpolating 
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between two filters, such as between / and K to get the iJ— band number counts, the 
results are less model dependent when the same model fits the counts in both filters in 
which they have been measured. I run the model in the old and new filters, and determine 
the predicted median color as a function of magnitude in the first filter. I subtract the 
median color from the data magnitude, and determine the ratio of the measured counts to 
the predicted counts. I then multiply the model prediction for counts in the second filter by 
this ratio, and this is the prediction for the number counts in the second filter based upon 
the measured counts in the first filter. The error bars in log{n) are the same. 

I will begin with a relatively simple example, converting the WFPC2 U^qq number 
counts in the HDF to the ground-based Uj band, so that they can be plotted in the same 
figure. The top panel in Figure |^ is the ground-based Uj counts plotted with a model 
for the Uj counts (solid line), and the f/300 counts (dashed line.) To construct a model 
that fits the Uj counts, 1 took a qq = 0.02 universe, included the effects of evolution and 
dust, and then set the mix of galaxy types to match the type- dependent results of the 
MDS (Glazebrook et al. 1995; Driver et al. 1995) and the HDF (Abraham et al. 1996) 
in the /—band. This resulted in a late-type dominated galaxy mix with a higher overall 
normalization than the luminosity function measurements of Loveday et al. (1992). Over 
the intermediate magnitudes, the model predicts different slopes for the two filters, as well 
as an offset of about 0.5mag. The next panel in Figure |^ is the HDF U300 number counts. 
I have plotted three versions of the counts, including the counts measured in isophotal 
and total magnitudes by Williams et al. (1996), and the counts measured by Metcalfe et 
al. (1996), who used aperture photometry in the f/300 image at the positions of objects 
identified on the image. The latter measurement goes deeper than the former, but 
is potentially subject to a selection effect against objects with extremely blue f/300 — -B450 
colors. The next panel in Figure |^ is the model prediction for the median Uj — f/300 as a 
function of Uj magnitude. Finally, the f/300 counts have been transformed into f/j, and all 
the counts are plotted together. 

To trace a single data point as an example, the Williams et al. (1996) measurement 
of the counts at Usqqab = 25.75 is login) = 4:.8Smag~^ deg'"^. Converting from the AB 
magnitude system to the Vega magnitude system is an offset of —1.42 for the f/300 filter, 
which puts this point at f/300 = 24.33. The median f/300 — Uj color at this point is —0.71, 
so the equivalent point is Uj = 25.04. The model predicts that at f/300 = 24.33 the count 
will be log{n) = 4.59mag~^ deg""^, which is a 0.29 less than the measurement. The model 
predicts that the count at Uj = 25.04 is log{n) = 4.71mag~^ deg'"^. Thus the converted 
data point is log{n) = 5.00mag~^ deg~^ at Uj = 25.04. 

I model the number counts in the NICMOS ifi.o (F160W) filter in Figure ^ I have 
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used the open universe model, including evolution and dust which gives a reasonable fit 
to the /—band number counts plotted in Figure ^. This model also fits the i^— band 
number counts, and I have plotted the J— band counts in the top panel of the figure and the 
i^'— band counts in the bottom panel. The i^'— band counts have been transformed to the 
ifi.6 filter by the method discussed above, and the results are plotted in the middle panel. 
When background limited, NICMOS is expected to obtain a 5a detection of point sources 
with Hi Q = 25.0 in a one hour exposure. Figure ^ demonstrates that multi-orbit exposures 
with NICMOS will reach fainter magnitudes than current ground-based near-IR surveys. 



3.2. Extrapolation into the UV and Thermal IR 

The method described above can be used to extrapolate current observations into the 
UV and thermal IR, so long as appropriate choices of SEDs and evolutionary models are 
made. The GISSEL96 model extends into the UV and IR, although as discussed above it 
does not include the effects of dust internal to the galaxies. Dust absorbs the UV flux and 
re-radiates the energy in the mid- and far-IR. While the simple dust model described in 



Section p.6| works reasonably well at modeling the effects on the UV flux from galaxies, 
the mid- and far-IR is much more complicated. Galaxies have dust with a range of 
different temperatures, sizes and compositions, all of which affect the IR SED. In addition, 
starbursts due to merging events can dominate surveys conducted in the far-IR, and must 
be considered a separate population. The application of the techniques of this paper to a 
model of galaxy evolution in the far-IR will be considered in another paper. I consider 
here extrapolation of current observations into the UV, and into the mid-IR. While the 
mid-IR is also affected by some of the considerations discussed above, emission from faint 
galaxies at high redshift in the mid-IR comes from the rest-frame near-IR, a wavelength 
region well modeled by the GISSEL96 models. Thus the model is more appropriate for 
predicting the faint observations planned by SIRTF and NGST, than for predicting the 
brighter observations of ISO and WIRE. 

Figure |^ shows normalized tracings of the STIS Near-UV and Far-UV MAMA 
detectors in their clear modes, compared to the WFPC2 f/300, ground-based Uj and the 
WFPC2 i?45o filters. In Figure |T0| I have used ground-based f/— band observations and the 
WFPC2 U30Q observations of the HDF to predict the number counts to be seen with STIS. 
Planned observations of an area within the Hubble Deep field using the STIS UV detectors 
are expected to reach fainter than 28mag in NUV, and 25mag in FUV; the NUV exposure 
will reach fainter than the U300 image and will strongly constrain galaxy evolution in the 
UV. 
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I plot an extrapolation of the i^— band number counts into the mid-IR in Figure 0. 
Deep exposures at 3.5/im and 8/im are planned as part of the Legacy Science Program for 
the Space Infrared Telescope Facility (SIRTF) (Clemens, Greenhouse & Thronson 1996). 
Observations of the HDF have been made at 6.7/im with the Infrared Space Observatory 
(ISO) (Oliver et al. 1997), and the results are plotted in the middle panel of the figure. 
The planned observations by SIRTF will reach slightly deeper at 3.5/im than ground-based 
i^'— band observations, although the planned observations at longer wavelengths will not. 
It should be noted, however, that an extrapolation in wavelength like I have made in this 
figure is mo del- dependent, and actual observations will constrain the models. In addition, 
of course, reaching fainter magnitudes or higher redshifts is not the only observational 
goal, and the mid-IR promises to play an important role in understanding the effects of 
dust and star-formation at low to intermediate redshifts. At the highest redshifts and 
faintest magnitude, the Next Generation Space Telescope, a 6 to 8 meter passively cooled 
instrument tentatively scheduled for launch around 2007, (Stockman 1997), will resolve and 
study the rest-frame UV, optical and near-IR at very high redshifts. 

3.3. Color Distributions 

The color distribution of the galaxies in a survey holds more information about the 
state of evolution than the number counts. In the rest-frame, younger galaxies generally 
are bluer than older galaxies, but the K-correction can counter this effect. Comparing 
galaxies to evolutionary tracks without considering the number with each color will give 
an indication of redshift and evolutionary state for individual galaxies, but does not give a 
clear indication of the state of the entire population. Therefore, the color distribution as a 
function of magnitude must be modelled in order to interpret the data. 

Figure |T2| is a plot of the number of galaxies as a function oi I — K color for three 
i^— magnitude bins. The bright data are taken from Gardner et al. (in preparation) and the 
fainter data are taken from Gardner (1995). The models plotted in Gardner (1995) were 
based upon the original Yoshii & Takahara (1988) models which do not include the effects 
of dust. The primary effect of dust in this figure is to block the bright ultraviolet flux 
from rapidly starforming galaxies at high redshift which come to dominate the population 
at faint levels. Galaxies appear red in the models because of their high redshift. Model 
predictions of the colors of galaxies are often not smooth functions because of the discrete 
nature of the models, so in this figure I have convolved the output with a 0.2mag Gaussian 
filter to mimic the effects of photometric noise. While it is always possible to run the model 
with a higher resolution in redshift or magnitude, models usually use discrete types. In the 
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model described here, I have removed this restriction, and allowed interpolation between 
the different integral types. To fit the data I plot an open universe model with evolution, 
dust and dwarf galaxies. This figure is an indication of what this model can do, but the 
interpretation of color distributions is more subject to uncertainties in the free parameters 
than are number counts. 



3.4. Luminosity Function Conversion 

The three parameters of the Schechter luminosity function (LF) are highly correlated, 
and any comparison of luminosity functions measured in different surveys must take this 
correlation into account. The maximum likelihood method for determining a parametric 
LF from a field galaxy redshift survey without binning the data was originated by Sandage, 
Tammann & Yahil (1979). This method has the advantage over previous methods, (see 
Felten 1977 for a review), that it is relatively insensitive to the effects of clustering, using 
the assumption that the LF is independent of environment. The maximum likelihood 
method, however, is subject to a normalization constant that must be determined in 
another way, so it returns M* and a, and the error ellipse for those two parameters, while 
0* must be determined independently. It is possible to determine 0* by normalizing a model 
of the galaxy counts based upon a measured M* and a to observed bright galaxy counts 
(Mobasher, Sharpies & Ellis 1993; Gardner et al. 1997). Unless model counts are run for 
every point within the error ellipse in M* and a, however, this method does not return the 
correlation between the errors in all three parameters. 

To compare LFs determined in different filters, it is best to compare M* and a with 
a translation between the filters, and then to compare 0* by comparing the bright galaxy 
counts directly with the galaxy count model. The galaxy count model presented in this 
paper converts LFs between filters using a method developed by Yoshii & Takahara (1988). 
The type-independent or type-dependent values of M* are converted by adding the model 
rest-frame color, and a is assumed to be the same for each filter. This will sometimes 
have the effect of converting a type-independent LF into a type-dependent LF, when the 
rest-frame colors are type-dependent. To compare the error ellipses directly would require 
doing this procedure for every point within the error ellipse. However, since the sum of 
several Schechter LFs (i.e. one for each galaxy type used in the model) is not likely to equal 
a LF which is easily parameterizable as a Schechter LF, this is not a straight-forward task. 
It is much easier to compare LFs graphically. 

The top panel in Figure |13| is a comparison of the 6 j— band luminosity function 
measured by Loveday et al. (1992), with the band LF measured by Gardner et al. 
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(1997), converted to the 6j— band using this technique. The excess of bright galaxies in this 
figure is due to the assumption that is independent of type, causing the bluest galaxies 
to appear too bright in bj. This assumption has little effect on model predictions of the 
near-infrared properties of galaxies due to the low number of blue galaxies, but is magnified 
by the conversion of the LF to the blue filter. It is an indication of the uncertainties 
introduced into this type of model by the extrapolation over a factor of five in wavelength. 
The middle panel contains the measured i?— band LF of Lin et al. (1996), and the converted 
6j— band and band LFs. The bottom panel contains the i^— band LF and the 6j— band 
LF converted into K. The assumption that M^^ is independent of type is more accurate 
than in K. 



4. Summary 

Galaxy counts are one of the classical cosmological tests, but their interpretation 
remain difficult. Model predictions of the counts, and the color and redshift distributions 
of galaxies are subject to uncertainties in the spectral energy distributions and evolution of 
galaxies, and in free parameters specifying the luminosity functions, cosmological geometry, 
the number and distribution of galaxy types, and the effects of dust and merging. However, 
there is a growing body of observational data that constrains these models. In this paper, I 
have presented a model in which the free parameters can be easily adjusted, and which can 
be extended to include new parameters. The greatest power of this model is in converting 
observational data taken in one filter into another filter to compare with other data, and 
extrapolating to longer or shorter wavelengths. 

The software, parameter files and data presented in this paper are available at 
http: / / hires.gsfc.nasa.gov/ ~gardner /ncnioB| . 

I wish to acknowledge funding by the Space Telescope Imaging Spectrograph 
Investigation Definition Team through the National Optical Astronomical Observatories, 
and by the Goddard Space Flight Center. 
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Fig. 1. — The dependence of the fiducial /814 model on cosmological geometry and evolution. 

In each panel I plot the total number counts of galaxies (solid line), and the number counts 
of each of 5 types: E/SO (short dashed), Sab (long dashed), Sbc (dash-dot), Scd (short 
dash-long dash) and Irr (dotted). Also plotted are a compilation of the /—band data (open 
points) including the HDF (filled points). The panels on the left include the effects of passive 
luminosity evolution, the panels on the right are no-evolution, or pure K-correction models. 
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Fig. 2. — The dependence of the models on the redshift of galaxy formation. The data and 
the Zform = 15 models are the same as the evolution models in Figure |1|, and I have included 
models with Zform = 5 and Zform = 2.5. The effect of a lower redshift of galaxy formation is 
to introduce a turnover in the counts. The faint end of the counts is then dominated by the 
slope of the faint end of the luminosity function. 
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Fig. 3. — The effects on the models of internal absorption by dust. The models plotted are 
the same as the evolution models in Figure but now the effects of internal absorption by 
dust are included. The dust absorption reduces the strong UV flux that is predicted for 
forming galaxies by the population synthesis models, and eliminates the strong up-turn in 
the counts seen in the dust-free models at intermediate magnitudes. 
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Fig. 4. — The effects of merging on the models. The models include self-similar merging 
with conservation of luminosity density. Merging has the effect of lowering the luminosity of 
individual galaxies while increasing the numbers at fainter magnitudes. The high observed 
counts at faint magnitudes do not rule out the go = 0.5, \q = 0.0 model if merging is 
included. The top panel uses the formulation of Rocca-Volmerange & Guiderdoni (1990), 
while the bottom panel uses the formulation of Broadhurst et al. (1992). 
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Fig. 5. — The effects of an extra population of dwarf galaxies on the models. In each panel, 
the solid line is the total number of galaxies, the dotted line is the number of dwarf galaxies, 
and the dashed line is the rest of the galaxies. A population of dwarf galaxies with a steep 
faint end slope has very little effect on the bright number counts, but comes to dominate at 
faint magnitudes. This is because the large galaxies with faint apparent magnitudes are at 
high redshift, where there is little volume and thus low numbers, while the dwarf galaxies 
are at low redshift where there is greater volume. 
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Fig. 6. — The combined effects of dwarf galaxies plus dust, and merging plus dust on the 
models. 
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Fig. 7. — The translation of number counts from the ground-based U band to the U300 filter. 
The top panel shows the ground-based [/—band counts and a model which fits them (solid 
line), along with the model prediction for the C/300 counts. The next panel shows the HDF 
C/300 counts along with the model. The next panel shows the models median Uj — C/300 color, 
which is then used to translate the C/300 counts into the Uj filter, as explained in the text. 
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Fig. 8. — The interpolation of number counts between the / and K bands, into the NICMOS 
i?i.6 filter. NICMOS is expected to reach S/N=5 for Hi q ^ 25 for a point source in a 1 hour 
exposure. This is well above the model prediction for 1 object per field of view with camera 
2. 
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Fig. 9. — Normalized filter tracings of various UV filter-detector combinations, including 
the STIS Far UV and Near UV MAMA detectors in clear mode, the WFPC2 U300 and ^450 
filters, and the ground-based U filter. 




Fig. 10. — A prediction of the number counts to be seen with STIS with the MAMA detectors 
in clear imaging mode. 
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Fig. 11. — A prediction of galaxy counts in the mid-IR. In the middle panel I plot the 
integral counts, and compare to the ISO observations on the Hubble Deep Field (Oliver 
et al. 1997). These 6.7 /im counts are only marginally inconsistent with the population of 
normal galaxies in my model, although the effects of dusty starburst galaxies, active galactic 
nuclei, and thermal re-radiation by dust are expected to become increasingly important in 
the mid-IR. Imaging with SIRTF is expected to become confusion limited at the sub-//Jy 
level, while NGST is expected to achieve nJy photometry. 
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Fig. 12. — A plot of the number of galaxies as a function oi I ^ K color. The model has 
been smoothed with a 0.2mag Gaussian to mimic the effects of photometric noise. The dwarf 
galaxy model was used to increase the number of blue galaxies and dust was used to decrease 
the number of red galaxies. 
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Fig. 13. — Comparison of luminosity functions determined in the 6j— band, the i?— band, 

and the i^— band. The top panel contains the &j— band LF of Loveday et al. (1992) with the 
fC— band LF of Gardner et al. (1997) converted to the bj band using the methods described in 
the text. The middle panel contains the i?— band LF from Lin et al. (1996) with the 6j— band 
and fC— band LFs converted to the R band. The bottom panel contains the iT— band LF of 
Gardner et al. (1997), along with the 6j— band LF converted to the K band. 



